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Abstract 

We apply the quantum renormalization group to construet a holographic dual for the U(N) 
vector rnodel for complex bosons defined on a lattice. The bulk geometry becomes dynamical as the 
hopping amplitudes which determine connectivity of space are promoted to quantum variables. In 
the large N limit, the full bulk equations of motion for the dynamical hopping fields are numerically 
solved for finite systems. Fronr finite size scaling, we show that different phases exhibit distinet 
geometric features in the bulk. In the insulating phase, the space gets fragmented into isolated 
islands deep inside the bulk, exhibiting ultra-locality. In the superfluid phase, the bulk exhibits a 
horizon beyond which the geometry becomes non-local. Right at the horizon, the hopping fields 
decay with a universal power-law in coordinate distance between sites, while they decay in slower 
power-laws with continuously varying exponents inside the horizon. At the critical point, the bulk 
exhibits a local geometry whose characteristic length scale diverges asymptotically in the IR limit. 
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I. INTRODUCTION 


According to the AdS/CFT correspondence[l-3], quantum field theories are dual to grav- 
itational theories in a spacetime with one higher dimension. The extra dimension in the bulk 
can be interpreted as a length scale in the renormalization group (RG) how[4-8]. However, 
the flow generated along the extra dimension in holography is different frorn the conven- 
tional RG flow because the bulk theories are in general quantum theories. In order to 
rnake the connection between RG and holography precise, one has to introduce quantum 
RG (QRG)[9-11]. According to QRG, general ZTdimensional quantum held theories are 
equivalent to (D + l)-dimensional theories that include quantum gravity, provided that the 
quantum held theories can be regularized covariantly in curved spacetime[10, 11]. 

Conventional renormalization group describes the flow of the couplings (sources) gener¬ 
ated by coarse graining dehned in the space of all sources allowed by symmetry. In QRG, 
only a subset of the sources is included. The sources in the subset are promoted to quantum 
mechanical operators, and a quantum Hamiltonian governs the evolution of the dynamical 
sources under scale transformations. Quantum huctuations in the RG flow precisely capture 
the effect of other couplings that are not explicitly included. In the bulk, this amounts to 
the fact that composite operators are generated after lowering the position of the UV bound- 
ary by integrating out bulk degrees of freedom[7, 8]. Intuitively, sources becorne quantum 
mechanical because high-energy modes act as huctuating sources from the point of view of 
low-energy modes [9-12] . 

Despite the formal mapping between general quantum held theories and gravitational 
theories, classica! and local geometries are expected to arise only for a special set of quantum 
held theories [13-15]. In the presence of a large number of degrees of freedom, the quantum 
huctuations of the RG how becorne weak, which allows one to use a classica! description. 
However, it is much more non-trivial to have locality in the bulk. Generally, there exist 
inhnitely many non-local operators, and locality arises only when the sources for the non- 
local operators are suppressed. In QRG, the condition for the emergence of local geometry 
translates into stringent constraints on the beta functions of the quantum held theory[16, 17]. 
From the constraints, one can try to hnd quantum held theories that exhibit a local geometry 
in the bulk. On the other hand, one can use varying degrees of locality that emerge in the 
bulk as a diagnostic that differentiates one phase from another. The bulk locality serves as 
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a useful order parameter for characterizing phases of matter using geometry, which is one of 
the primary goals of applying holography to condensed matter systems[18-21]. 

In this paper, we apply QRG to a three-dimensional U (N) vector model regularized on a 
Euclidean lattice, which can be viewed as a quantum Bose-Hubbard model of N components 
in two space dimensions and an imaginary time. It supports an interacting field theory at the 
critica! point between the insulating (gapped) phase and the superfluid (long-range ordered) 
phase. The full bulk equations of rnotion derived from QRG are solved in the large N limit 
on finite size lattices. From a finite size scaling analysis, we show that different phases indeed 
exhibit different degrees of bulk (non-) locality in the thermodynamic limit. 

Here is the outline of the paper. In Sec. II, we start witli the U(N) vector model for N 
complex bosons dehned on a D-dimensional lattice. In this model, the general single-trace 
operators can be written as 0* a 0j a , where (j) ia is a complex boson with flavour a dehned 

at site i. The single-trace operators involve only one flavour contraction and are bi-local in 
space[22-24]. The set of single-trace operators has a special status in that ali other singlet 
operators can be constructed as composites of single-trace operators. They describe the 
kinetic terni that connects pairs of sites through hopping. In the vector model, all single-trace 
operators are quadratic, and it is not enough to include only single-trace operators to describe 
the phase transition from the insulating (gapped) phase to the superfluid (symmetry-broken) 
phase. Therefore, we also include a quartic double-trace operator. We treat the single-trace 
hopping ternis as deformations to the on-site action that includes the single-trace mass terni 
and the double-trace interaction terni. I 11 the conventional RG, a series of coarse graining 
generates an infinite tower of multi-trace operators [25]. In QRG, one only keeps track 
of single-trace deformations at the expense of promoting the sources for the single-trace 
operators to quantum variables. We derive a (D + l)-dimensional bulk action whose degrees 
of freedom are scale dependent bi-local hopping helds (■ Uj{z )) and their conjugate variables 
(t*j(z)), where z is the logarithmic lengtli scale. The bulk theory describes the quantum 
mechanical RG flow in the space of single-trace operators. In the large N limit, the path 
integral in the bulk can be replaced by saddle point equations whose solution determines 
the classical geometry in the bulk. In Sec. III A, an analytic solution to the bulk equations 
of motion is found in the deep insulating phase by treating hoppings as small perturbations 
compared to the on-site ternis. Away from the deep insulating phase, we resort to numerical 
Solutions. Sec. III B and C outline the numerical scheme. I 11 Sec. IV, the numerical 
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Solutions obtained for three dimensional lattices with linear sizes 3 < L < 13 are presented. 
Tlien we extract the behaviour of the solution in the thermodynamic lirnit frorn finite size 
scaling. 



FIG. 1: A schematic holographic phase diagram in the thermodynamic limit. A rh 2 is the mass 
which tunes the System from the superfluid phase to the insulating phase. z is the extra dimension 
in the bulk, which corresponds to a logarithmic length scale in the renormalization group. In the 
insulating phase, the space gets fragmented beyond the fragmentation scale, zp. In the superfluid 
phase, non-locality emerges beyond the horizon scale, zh • At the critical point, the bulk exhibits 
a local geometry. 

Fig. 1, which is the key resuit of the paper, summarizes the numerical solution in the 
thermodynamic limit, which is suggested from the finite size scaling analysis. Here A m 2 
represents the on-site mass term, which tunes the systern away from the critical point either 
to the insulating phase (Arh 2 > 0) or to the superfluid phase (Arh 2 < 0) depending on its 
sign. Although only the nearest neighbor hoppings are turned on at the UV boundary at 
z = 0, further neighbor hoppings are generated inside the bulk with z > 0. In the insulating 
phase, the hopping fields Uj(z) decay exponentially not only in |i — j | but also in z. As a 
resuit, the lattice is fragmented into decoupled sites in the large z limit. While the scale Zp 
beyond which fragmentation occurs remains small in the deep insulating phase, it diverges 
as the critical point is approached in the thermodynamic limit. In the superfluid phase, on 
the other hand, a new scale emerges. Instead of fragmentation, the systern loses locality 
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beyond a critica! scale zh, which is called the horizon. At the horizon, the hopping fields 
decay in a universal power-law in \i — j\. Inside the horizon with z > Zh, the hopping fields 
decay with a slower power-law with continuously varying exponents. The horizon scale zh 
diverges as the critical point is approached frorn the superfluid side. At the critical point, the 
bulk exhibits a local geometry with a characteristic length scale that remains finite for finite 
z and diverges only in the IR limit. In summary, the insulating phase, the superfluid phase 
and the critical point exhibit ultra-local, non-local and local bulk geometries, respectively. 
In this sense, locality serves as a holographic order parameter that characterizes different 
phases of matter. The goal of Sec. IV is to justify the holographic phase diagram in Fig. 1. 
For related works on the holographic description of the vector models, see Refs. [22-29]. 


II. HOLOGRAPHIC ACTION FOR VECTOR MODEL ON A LATTICE 


We consider a bosonic vector modcl in H-diinensional Euclidean space, 


S 


d D x 


\X7(p\ 2 + m 2 \(f>\ 2 + -^r(|0| 2 ) 2 


( 1 ) 


where <f> refers to the complex boson field with N flavours, \(f>\ 2 = <f>* ■ <f>, m is the mass 
of the bosons, and A is the quartic coupling. We regularize the continuum theory on a 
H-dimensional hypercubic lattice, 

§ = S 0 + S \, (2) 


where 

50 = m 2 Y^ (<f>* ■ <t>i) + , ( 3 ) 

i i 

51 = - 4' (<t>i ■ 4>j) + Yl (^i ‘ Qj) (<t>*p • <t>q) • ( 4 ) 

ij ijpq 

Here i,j,p,q run over the ZVdimensional Euclidean lattice for the N component bosonic 
fields (fi a with a = 1, • • • , N. So represents on-site ternis, and Si includes the hopping terni 
and general quartic interactions. To guarantee the action is real, we impose f t j = tf and 
Jijpq — Jjiqp- Here we consider the case where f tJ is real. We could make f tJ complex by 
turning 011 a background gauge field. The background gauge field in general breaks the 
time-reversal symmetry, if we go to Minkowski space by Wick-rotating one of the Euclidean 
directions into real time. For general tij and Jij pq , the model has the U(N ) symmetry. 
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The symmetry is enhanced to 0(2N) when t t j's are real and Jij pq oc SijS pq . The action is 
proportional to N in the large N linxit with fixed m 2 , A, t^ and Jij pq - 

There is a redundancy in the parameters of the action, and the action is invariant under 


77i —y Tn T A, tij —y tij -t- ASij, 


A y A 13. Jijpq ^ Jijpq BdijSipSiq. 


(5) 


By shifting tu and J lul , So coulcl liave been entirely absorbed into S\. Here So is explicitly 
singled out because we will use So as ‘the reference theory’ and treat S± as a deformation. 
In other words, the renormalization group (RG) flow will be dehned in terms of the flow 
of the parameters in S± which is dehned with respect to the fixed So with m 2 , A > 0. 
It is emphasized that one can describe not only the insulating (gapped) phase but also 
the superhuid (symmetry broken) phase because t l j can be arbitrarily large to support a 
Mexican hat potentia! for (fi. The freedom to choose different m 2 and A amounts to choosing 
different RG schemes, which does not affect physica! observables. 

We derive a holographic theory for the lattice action in Eq. (2) using the QRG scheme. 
We start by writing the partition function: 


Z = / T)<fiT)(fi*e 


* —s 


where the measure is given by T)<fi = n=i u!Li d(j) ia , V<fi* = nL nf =1 dfi* a and V is the 
number of sites in the lattice. The hrst step is to remove multi-trace operators in Sf by 
promoting the sources for the single-trace operators into dynamical variables. Using the 
representation of the Dirac-Delta function[9] 


N 


f(A ) = — dt dt* e -^ m *- A) f(Nt*), 


7T 


( 6 ) 


we introduce a pair of complex conjugate link helds, t\j and t*j for every ordered pair of 
sites i, j to rewrite the partition function as 

Z= [ D0H0*Rt (o) Rr (o) e- s(o) , 


where Vt^Vt*^ = U uj dtff dt*f 


with i,j running over ali sites and 


§(°) = So + NS 


uv 


AO) ,*(0) 
°ij ’ L ij 


XX 0) (^i ' & 3 


(7) 
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with 

= - E kC +E 4« C C '+E CC- (8) 

y ijPQ ij 

Here and henceforth, multiplicative constants in the partition function are ignored. In the 

new action, only single-trace deformations are present. Although the donblc-trace operator 
in Eo could have been removed as wcll, we need to keep the on-site donble-trace operator to 
make sure that the patii integra! for <p is well defined for any values of fluctuating sources . 
The presence of the double-trace operator is also important because it rnakes the hopping 
ficlds genuinely dynamical variables. 


~>uv 


ao) t m] 


y 


y 




(a) 


(b) 



(c) 


FIG. 2: (a) A lattice of four sites where every site is coupled with all other sites through dynamical 
hopping fields. (b) Upon permutation among a subset of sites represented by solid (red) arrows, 
the action remains invariant if the dynamical hopping fields are also permuted as is denoted by 
dashed (blue) arrows. (c) An example of a permutation that doesn’t preserve the local coordinate 
volume element. The blue lattice sites remain fixed, while the red lattice sites move according 
to the denoted arrows. Under the transformation, some rectangles with coordinate volume 1 are 
rnapped to regions with coordinate volume 1.5 or 0.5, as is indicated with shaded regions in the 
figure. 


Note that <fi is coupled to dynamical sources tf^ which connect all possible sites. There- 
fore the new action is defined on a globally connected network rather than a fixed lattice. 
The globally coupled network with dynamical hoppings does not have a fixed lattice struc¬ 
ture, and i,j should be regarded as indices for ‘events’ rather than fixed coordinates. Since 
there is no fixed lattice, there is a freedom to relabel each event differently, which corre- 
sponds to a discrete local coordinate transformation. So we perform a discrete coordinate 
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transformation, 


<t>i ^: + n (0) ’ (9) 

where site i is displaced by a site-dependent D-dimensional shift vector, N^[10, 11, 26]. 
Since the full set of events rnust map to itself, the transformation shonld form cyclic per- 
mntations of the events. This is illustrated in Figs. 2(a),(b). After the local coordinate 
transformation, Sq is manifestly invariant and the hopping term is transformed to 


- ' <t>j) ( 10 ) 

ij ij 

where the shifted hopping held is given by 

/(o) = Ao) 

V i-Ni 0 ) ,j-N (0) ' 

* 3 

The point is that the change introduced by the local coordinate transformation can be 
compensated by a transformation of the hopping helds because tf-' is dynamical. Therefore, 
the dynamical hopping helds play the role of a dynamical metric in the continuum limit, 
which guarantees invariance under a set of .D-dimensional local coordinate transformations in 
the continuum limit. The set of permutations generated by the shift includes D-dimensional 
local coordinate transformations, which don’t necessarily preserve the coordinate volume 
locally. An example is given in Fig. 2(c). The fact that the hopping helds play the role 
of a metric is expected because the physica! distance between two sites is determined from 
the strength of the hopping between them: the larger the hopping between two sites, the 
shorter the physical distance. 

Now we generate RG how by coarse graining the System in real space. The original held 
cj) is split up into liigli- and low- energy helds, where the low-energy held has a mass slightly 
larger than m 2 . The missing huctuations are carried away by the very massive high-energy 
helds. Then, the high-energy modes are integrated out to generate quantum corrections for 
the low-energy modes[30, 31]. For this, we introduce an auxiliary held $ with an arbitrary 
mass /i, 



Z= I D</)D(/)*I)$D$*D^ (0) DD (0) e“ s(0), , 


where 

S (0) ' = NSw 


Ao) Mo) 

L ij ■> L ij 


4 0) (^< ■ H 2 ■ fc) + +4 ■ ^) 2 • 


N 




10 




Now we rotate the physica! field and the auxiliary field into high- and low-energy fields, 

0i — 0; + 4>i 

= Aify + Bityii ( 12 ) 

with Ai = Bi = — ^ and /I,; = , a\ ' is the ‘lapse’, which is a site-dependent 

M e i dz -l 

constant that Controls the local speed of coarse graining[10, 11, 32]. Here < fi' is the low- 
energy mode with mass m 2 e 2oii )dz , which is larger than m 2 by an inhnitesimal amonnt 
0(dz), and 0 is the high-energy mode with a large mass /lj 2 e 2< 0 )dz . After rescaling the 
fields 0 ' — > e~ a i ]dz 0(, 0j —> e _ "i )<iz 4>i and renaming 0' —>■ 0,., we rewrite the partition 
fnnction as 


Z= H0H0*I)0H0*2)t (o) T)r (o) e- s(o),, 


where 


§(°)" = N§uv + Y p 2 (0* ■ 0O + /h 

i 

( 0 * + 0 t) • (p + 0 . 


2 / j * 


4 >i 


- Vt (0) e -(“‘ 1)+a i 1) ) ,iz 


A_ 

IV 


E' 


—4q^ 1 ^2: 


0* + 0*) • 0* + <t>i 


1 2 


(13) 


Integrating out the high energy mode 0, we obtain the renormalized action for the low 
energy mode, 

1 /(o) 


§ (0) " =iV Sj/y [4 0) ’4' 0) J + 2AhfeJ]cq (1) | t 


+ 2ch ctj 


(i) f 2A (l + 


rrr 


( 0 * • 00 


4A 2 


m 2 N 2 


(0* • 0t) :E 


+ 2 dz Y <{ (4 O) (0i ' 0i) + 0f(0j ■ 00) (0*' 00 

+ 2 * E «f’ { -i E ‘2W« ■« - 4 <«' + E ^ w ■ *) til (*? •*■) 


rrr 


jk 


Y f T ( 0 i ■ 0 j) + m2 XI ^ ’ 


A 


AT 


(14) 


*7 


where quantum corrections are kept only to the linear order 0(dz) (see Appendix A for 
details). In the renormalized action, multi-trace operators have been generated for the low 
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energy field. Under the conventional (classica!) Wilsonian RG procedure, one would apply 
the same coarse graining procedure to Eq. (14) to generate even higher-trace deformations. 
In QRG, we remove the multi-trace deformations before repeating the coarse graining. Tliis 
way, one can keep only single-trace deformations at each step of coarse graining by making 
the sources for the single-trace deformations fluctuating. Just as before, we introduce new 
link fields t^] and £*,■ to remove the multi-trace deformations, 

Z = I V4> Txfi* T>t {0) RU (0) Dt (1) Vt* w e" s(1) , 

where 


s (1) =" Suv + N - *T ) d 1 ’ 




m 2 Zj bj L kj 


- 2A U 


*(i) 


E 


/(o) ,*(i) , /(o) ,*( 1 ) 

L ij L ij i L ji L ji 


O # • fa )+ m2 + ■ ^) 2 • 




A 


iV 


v 


(15) 


Now we can repeat the coarse graining procedure as before. The coordinates of the cfifis 
can be shifted by a new set of discrete translations, N^b Then (fi is split into low-energy and 
high-energy modes where the low energy mode at site i has a mass m 2 e ai dz . After rescaling 
the fields, we integrating out the liigh energy modes to generate quantum corrections which 
include multi-trace operators for the low energy fields. We introduce another set of auxiliary 
fields and t*^ to remove the multi-trace operators. From here the pattern is ciear. After 
doing this T times, the partition function is written as 

Z= f VcfiVcfi* JjRf (/) RU w e- Sr , 

^ 1=0 
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where 


Sr —N S 


uv 


'C.Cl + N E 


Z=1 




2 *E“< 


L 

+ 


1 /O-i) 2A (l + ,* m 4A 2 


m- 


mi ) _ (EE Mo 

“ m 2 V * 


2A /(i-i).*(i) , /0-iR*0)w*0) _ 1 WQ-i)/Q-iy 
+ m 2 Zj' bj tl i< L ji ' L ii m 2 Zj k bj 1 


—1)+*0) 
kj 


jk 


- 2A t. 


L *(0 


/O-i) ,*0) , /0-i) ,*0) 

L ij L ij ' L ji L ji 


4 r) (# • 0 i) + m2 S ^ (# ' &) 


(16) 


with tfp = t® (J) . (!) . Here i®, are the dynamical sources and thcir conjugate 

variables introduced at the Z-th step of coarse graining. The site-dependent lapse and shift 
{cxf^N^} are chosen independently at each step of coarse graining. Finally, we take the 
continuum limit in the RG direction by introducing z = Idz in the limit dz —> 0 with 
£* = Ydz hxed . In the limit where 0 is continuous, tf), tf- \ af \ N ; !/) become tij(z), t*-(z ), 
cti(z), Nj(z). The partition function is written in terms of the path integra! for the scale 
dependent hopping helds and their conjugate helds, 


2 j — r Dt r Dt* e~ N ^ uv+Sbulk+ ^ IR ^ 


(17) 


Here §uv , dehned in Eq. (8), is the UV boundary action for the dynamical source at z = 0. 
Sbuik is the bulk action, 

Uj{z) - tUz) 


§>Buik= / dzl YMz) d z tjj(z) + 




dz 


y2^i(z)\ — -t u (z )+ 

\ m 


2 _ 4A(1 + .8A 2 


m- 


<i(z)~ 4A (tl(z)Y- — (tl(z)y 


m- 


4A 


+ 534 ^) f ij( z ) +t 'ji( z ) f U z ^ + —2 + t 'ji( z ) t *ji( z ^ t *n( z 

777 / 


777 / 


Ehi( 2 )4w f tMi 


(18) 


jk 


with zh(^) = Zi-N i (z),j-N ■.,•(*) ( 2 )- SVb is an action, 


SjjiM**)] = - ]y ln 


T>cj)T>(j)*e~ 


-[- £ij Uj{z*)(<l>*-<l>j)+m 2 Ei(^i+0+W 


(19) 
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which is defined at the IR boundary z = z*. z* is the scale at which we stop the coarse 
graining procedure, which can be taken to be infinite. 



FIG. 3: Two hopping fields t^i and tij fuse to generate a longer range hopping between sites k 
and j. 

The new action is written in terms of the fields t t] {z) and t*j(z) which are bi-local within 
the original D-dimensional lattice and local in the 0 direction. If one interprets z as ‘time’, 
the action describes a System of bi-local hopping fields under a Hamiltonian evolution in 
imaginary time, tij and t*j correspond to annihilation and creation operators respectively[9], 
and the last term in Eq. (18) describes the process where two hoppings merge to become 
a longer-range hopping. This is illustrated in Fig. 3. In the path integra! language, ali 
RG paths are summed over for the single-trace operators. The weight for each RG path is 
determined by the bulk action, Ssuik- 

In the bulk, the lapse and shift vector parameterize local coorclinate transformations. The 
partition function does not depend on {eti(z), Nj(z)} because different choices of lapse and 
shift rnerely correspond to different renormalization group schemes [10, 11]. As a resuit, the 
local coorclinate transformation in the bulk is a gauge symmetry. Unlike in the continuum, 
the action contains terms that are non-linear in the shift. This is because the discrete shift 
is a finite transformation. When Uj(z) vary slowly in i,j, one can take a continuum limit 
along the H-dimensional directions by writing 

tij(z)->t(r 1 ,r 2 ,z), (20) 
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where ri = arj, r 2 = arj in the limit that the lattice spacing a is small. Then, the shift 
can be performed infinitesimally, Nj(^) —> n(r ,z)dz, and one can drop the terms that are 
non-linear in dz. In this case, the action is linear in n(r, z) and a(r ,z), which generate 
(D + l)-dimensional coordinate transformations in the bulk. Also, in the continuum, the 
bi-local helds can be represented by an infinite set of helds with arbitrarily large spin. Here 
we don’t take the continuum limit and proceed with the lattice action. 

In the absence of the quartic terni (A = 0), the theory has a larger symmetry generated 
by[28, 29] 


3 

l 3 » IJ 33 


./ ./ 
* 3 


where V %3 is a unitary matrix in the space of ali sites in the lattice. The enlarged symmetry 
is the source of the liigher spin gauge symmetry for the free theory with A = 0. In the 
interacting theory with A ^ 0, the symmetry is broken down to the local coordinate trans- 
formation. One might keep the liigher spin symmetry in the bulk action even when A 7 ^ 0 
by absorbing the quartic terni into the UV boundary action Suv- In that case, there is 110 
A in the bulk, and the quartic terni is implemented only through the alternative boundary 
condition determined by Suv [33] . However, this can be done only in the striet large N limit 
where the fluctuations of the dynamical hopping field induced by Suv are negligiblc. What 
is less satisfying in this description for the interacting theory is the fact that 1 /N corrections 
are singular because the path integral for <fi is ill-dehned without the quartic terni in Sq for 
fluctuating hopping helds. Therefore we choose to include the quartic terni in the bulk for 
general N, which breaks the liigher spin symmetry. 

Because the helds tij,t* 3 are singlets under U(N ) rotations, the action in Eq. (17) is 
proportional to N in the large N limit. Therefore, the partition function can be obtained by 
the sadclle point approximation in the large N limit. To that effect, ali we need to solve are 
the equations of motion, which we focus 011 from now 011 . Although t tJ and £*■ are complex 
conjugates of each other in Eq. (17), their saddle points conhgurations do not have to satisfy 
that condition. This is because the aniplitude and phase of tij can take complex values at 
the saddle point. Therefore, we treat and t*j as independent complex helds in the saddle 
point equation. 
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In the gauge with cti(z) = 1 and Nj(z) = 0, the equations of motion for the bulk fields 
are given by 


d z tij = -2 { ^t 13 - Sa 


rn 




' , 12A 2 

4A + — —t. 


m 


2 ** 


t ii + 


2A 5, 




m- 


'y ' (tiktik + tkitki) 


1 + m 2 + t jj) t% 3 2 tiktk i \ ’ 

J k J 

1 + ^2 (*»» + tjj) tij ~ rn2 X/ (^jk + f ) 


^ = 2 <| —^ + 

j i 


rrr 


m- 


( 21 ) 


subject to two boundary conditions 


iy( 0 )*i ii + 2 ^ 4 M t^( 0 ) = 0 , 


pq 


COO + = 0 , 




( 22 ) 

(23) 


which are derived by extremizing the boundary actions at z = 0 and z = z*, respectively. 
In the equations of motion, sub-leading ternis are dropped in the large N limit. From the 
IR boundary action, one can readily write the IR boundary condition as 



• 0 ,)>»,», 


(24) 


where < O >§ IR denotes the expectation value of O evaluated with respect to §jr at z = z*. 
This relationship must hold anywhere inside the bulk, since the RG process can be stopped 
at any point. Therefore, the on-shell value of t*j(z) coincides with the 2-point correlation 
functions at scale z, which are completely determined by Uj(z). 

In the end, one has to solve the hrst order differential equations for tij and t*j with one 
set of boundary conditions imposed at the UV boundary and the other imposed at the IR 
boundary. It is the IR boundary condition that imposes the constraints that the expectation 
values of operators (t*-) have to satisfy for a given set of sources (Uj). In particular, the 
vacuum expectation value has to satisfy the Ward identity associated with the discrete 
coordinate transformation, 


<0;-0i)i = W +N r0j+ Nj V’ < 25 ) 

where t is dehned in Eq. (11). Since this equation has to be satished at ali z, it becomes 
a dynamical constraint, t*j(t(z)) = t* +N .j +Jsl . (t'(z)), which is the discrete version of the 
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energy-momentum conservation imposed by N;(z). This condition is automatically implied 
by the IR boundary condition even in the fixed gauge with Nj(z) = 0. It is noted that t*- 
cannot be arbitrary at the IR bonndary, but it has to be the actual vacuum expectation 
value computed from the IR boundary action, Sjr. 


III. SOLUTIONS TO SADDLE POINT EQUATIONS 

In this section, we examine the solution to the equations of motion both analytically 
and numerically. We assume that translational symmetry is present, in which case t^ and 
t*j depend only on | i — j\ = |r* — rj\. Also, since we are maintaining the full SO(D) and 
.D-dimensional inversion symmetries at each scale z, all tij(z) and t*j(z) are real. We focus 
on the tliree dimensional cubic lattice of linear size L with periodic boundary conditions. 

To uniquely determine the Solutions we need to impose the boundary conditions in Eqs. 
(22) and (24). At the UV boundary, Eq. (22) imposes the Dirichlet boundary condition 
for tij( 0) = in the absence of J l]pq - To impose the IR boundary condition in Eq. (24) 
we compute §>m[tij(z*)] by introducing a Hubbard-Stratonovich field cr which satisfies the 
saddle-point equation in the large N limit (see Appendix B), 


a = 




V ^ m 2 + 2a — tJz*) ’ 
k&v ’ 


(26) 


where tk(z) = JT- e ' k ' (rj ri Hij(z). Then the IR boundary condition becomes 




m 2 + 2a — t k (z*) ’ 


(27) 


A. Analytic Solutions in the deep insulating phase 

For m 2 tij and m 4 A, the quadratic on-site action dominates. All correlation 
functions decay exponentially, resulting in the insulating phase. In this case, one can easily 
obtain an analytic solution to the equations of motion as a perturbation series in and 
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To order 1/m 6 , the Solutions are 


,* _ $ij_ 

ij m 2 
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tij — UjC 2z + 
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^ ' tik^kj 2X8, 
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)£«*, 


— (e 6z — 2e 4z + e 2z ) ^ ^ Ukikk'ik'k"ik"j • (29) 

fc /c' fc" 

Botli Uj(z) and t* 3 (z) decay exponentially as one moves towards the IR boundary. This is 
expected since bosons are localized at low energies. It is noted that t* 3 for i ^ j decays as 
e~ 2z witli increasing z. This feature holds generally even when t t3 is not small. This is due 
to the fact that the transformation in Eq. (12) followed by the rescaling by e~ dz implies 
4>i{z + dz) = e~ dz (pi(z ) np to the contribution from the anxiliary field introduced at each 
step of coarse graining which has zero correlation lcngth. This implies that 


< 4>*(z + dz) ■ 4>j(z + dz) >= e 2dz < <f>*{z) ■ <t>j{z) > 


(30) 


for i 7 ^ j. 

In the insnlating phase a large simplihcation can be made for the IR boundary condition. 
Since both Uj(z) and t* 3 {z) decay exponentially in z except for i = j, we can approximate 
them as Uj(z*) ~ and t*-(z*) ~ ^ for a sufficiently large z*. In this case, the IR 
boundary condition reduces to a single-site problem, 


t*-(z*) = _ — _ 

} m 2 + 2a — tu(z*) ’ 

where a is given by a = | (ta(z*) — m 2 + yj(tu(z*) — m 2 ) 2 + 8 A 
the on-site saddle point equation a = -g 0 x _, , , v 

772 “r Z(J tu yZ J 


(31) 

which is the solution of 


B. Numerical Solutions 

Since it is hard to solve the full equations of motion analytically in general, we resort to 
numerical Solutions for the superfluid phase and the critica! point. We solve the coupled 
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differentia! equations using a spectral method. The details of the numerical method can be 
found in Appendix C. Nnmerical Solutions are obtained for the cubic lattice with linear sizes 
3 < L < 13. We set A = 1, m 2 = 25, and A,- = 0 except for the on-site and nearest neighbor 
hoppings. The nearest neighbor hoppings are set to be 1, and the on-site ‘hopping’ tu is 
tuned to change the physical bare mass, 

rh 2 = m 2 — tu (32) 

to drive the phase transition from the insulating phase to the superfluid phase. The physical 
mass m 2 can be either positive or negative whereas the mass m 2 of the reference action in 
Eq. (3) is fixed to be positive. The reason for choosing a positive m 2 is to make sure that 
only short-distance modes are integrated out in the coarse graining procedure. 

Although there is no genuine symmetry breaking phase transition in finite Systems, fea- 
tures of the superfluid phase and the critical point can be stili identihed from the behaviors 
at intermediate length scales before the finite size effect takes over. Henceforth, we will keep 
using the terminology “superfluid” to refer to Systems in parameter regimes which become 
real symmetry broken states in the thermodynamic limit. Although only the nearest neigh¬ 
bor hopping is turned on at the UV boundary, ali further neighbor hoppings are generated 
inside the bulk. From now on, we will use the notation tn mn ) denote the specikc 

fields with r j — rj = ( l,m,n ) in units of the lattice spacing. For example, the three 

nearest neighbor hopping fields are represented by t( 0 i o,i )(^) 5 fyop.o)^), ^(i,o,o)(^)- 

C. Universal IR boundary condition in finite Systems 

In Figs. 4 and 5, we plot Uj(z) and t*j(z) for L = 3 with z* = 8. For each ( l,m,n ), 
we display multiple curves at 25 > m 2 > —15 that cover from the deep insulating phase 
to the deep superfluid phase across the critical point, which is at m 2 = 5.49454 in the 
thermodynamic limit. Since the curves look very similar away from critical point and they 
start to change rapidly only in its vicinity, we choose the spacing of m 2 to be smaller near 
the critical point and larger away from it. Here, we notice something interesting. No matter 
what the values of rh 2 are, ali Uj(z) and t* 3 (z) with i ^ j approach zero in the large z 
limit. In other words, the systern flows to the on-site problcm in the deep IR not only in 
the insulating phase but also in the deep superfluid state. The universal insulating (on-site) 
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*(0,0,0) ^(0,0,1) 



(c) (d) 


FIG. 4: The scale dependent dynamical hopping fields t^ >rnjn -\(z) for (l,m,n) = rj — rj = 

(0,0, 0), (0,0,1), (0,1,1), (1,1,1) for L = 3 and z* = 8. Different curves in each panel are at 
different values of m 2 = 25, 20, 15, 10, 7.5, 6.5, 6, 5.49454, 5, 4.5, 4, 3.5, 2.5, 0, -5, -10, -15. 
The hopping fields increase monotonically with decreasing rh 2 , as is denoted by the arrow in (b). 
The thick line is for rh 2 = 5.49454 which is the critical point in the thermodynamic limit. The 
curves below (above) the thick lines are in the insulating phase (superfluid phase). Signatures of 
the phase transition will be discussed in Sec. IV in detail. 

behavior in the deep IR limit is due to the finite size effect. For m 2 < m 2 , tij increases with 
increasing z for small z, which is consistent with the expected RG flow in the superfluid 
phase where the hopping terms become more important at larger lcngth scales. However, 
the RG lcngth scale eventually becomes greater than the System size. Then the finite size 
effect becomes important, and the systern “flows” back to the insulating behavior in the 
deep IR. This has to do with the fact that ali pliases are adiabatically connected to the 
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(c) (d) 


FIG. 5: The scale dependent correlation functions, ^ (z) for (l,m,n) = r j — rj = 

(0,0, 0), (0,0,1), (0,1,1), (1,1,1) for L = 3 and z* = 8. The same values of m 2 are chosen as 
in Fig. 4, and the magnitudes of the correlation functions increase monotonically with decreasing 

~ 9 

m. 

insulating phase in finite Systems. 

At first, the finite size effect appears to be detrimental to our goal of characterizing 
different phases of matter in terms of different IR geometries. However, one can stili study 
the superflnid phase and the critical point by focusing on the behavior of t VJ , t*- in the 
intermediate range of z before the finite size effect becomes dominant. Moreover, we can 
use the finite size effect to our advantage. The universal insulating behavior in the IR limit 
allows us to use the singlc-site IR boundary condition in Eq. (31), which is much simpler 
than the full IR boundary condition in Eq. (27). This has a significant implication. A 
practical difhculty in applying QRG to strongly coupled held theories is that one has to 
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impose the IR boundary condition dynamically. This is done either by extremizing the bulk 
action in the z* oo linxit or imposing Eq. (24) at a finite z*. If one uses the latter scheme, 
one has to know Sjr. Although this does not pose any problem for the present vector model 
which is exactly solvable in the large N lirn.it, this is in general a difficult task for strongly 
coupled theories (such as matrix models). However, if one can always impose the insulating 
boundary condition for finite Systems one only needs to know the solution to the on-site 
problem, which is much easier. This allows one to find Solutions for finite Systems without 
knowing the IR behavior of the infinite System a priori. Then, one can extract the behaviors 
in the thermodynamic limit through finite size scaling. This is the strategy we will employ 
in the rest of the paper. 

In order to guarantee that Systems with finite sizes flow to the insulating phase, one has 
to choose z* to be sufficiently large. The smaller m 2 is, and the larger L is, z* has to be larger 
to ensure that the system flows to the deep insulating state driven by the finite size effect at 
the IR boundary. This is because it takes longer RG ‘time’ before the finite size effect takes 
over in deeper superfluid phases and larger lattices. We observed that for all values of m 2 we 
took within the range 25 > m 2 > —15 and all lattice sizes 3 < L < 13, z* = 8 is sufficient. 
For larger L , it is expected that z* should grow logarithmically in L because z corresponds 
to a logarithmic length scale. Although computational cost increases with increasing z* } it 
is stili much cheaper to use the on-site boundary condition with a logarithmically larger z* 
than using the full boundary condition whose computational cost increases much faster. In 
Appendix D, we show that the Solutions obtained from the on-site boundary condition are 
indeed identical to the ones obtained with the full boundary condition. From now on, we 
will use the on-site boundary condition which allows us to obtain full numerical Solutions for 
larger lattices. The plots of the hopping fields for larger lattices can be found in Appendix 
D. 

IV. ANALYSIS 

In this section, we present the main results of the paper. We first establish the presence 
of the critical point that divides the insulating phase and the superfluid phase based on the 
correlation functions measured at the UV boundary. Then we examine how different phases 
exhibit distinet geometric features in the bulk. 
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A. UV Boundary 


^(4,4,4) 



FIG. 6 : The correlation function at the UV boundary t*-{z = 0) = -^ ((</>* • 4>j))uv between 
farthest sites (l,m,n) = r j — ri = (4,4,4) plotted as a function of rh 2 for L = 9. Although there 
is no true singularity due to finite size effect, the correlation function exhibits a rather sharp kink 
at the critical point, which is located at rh 2 ~ 5.5 in the thermodynamic limit. The straight line is 
Tjd&i ■ <t>j))uv = -hdS which is the expected value in the thermodynamic limit. 

The correlation functions, ■ <fij)) are given by t*j(z = 0), which are measured at the 

UV boundary. In the thermodynamic limit, lim T ._ >oo t* i+r (0) is the order parameter for the 
phase transition between the insulating (symmetric) phase and the superfluid (symmetry 
broken) phase. Although there is no real phase transition in finite lattices, there exists 
a rather sharp crossover which becomes the phase transition in the thermodynamic limit. 
To identify the ‘would-be’ phase transition in finite Systems, we plot i*-(0) for the largest 
possible \i — j\ as a function of rh 2 for L = 9 in Fig. 6. Although it is not strictly zero even 
in the insulating phase due to finite size effect, there is a sharp crossover around rh 2 « 5.5. 
The critical point in the thermodynamic limit is at rh 2 = 5.49454 (see Appendix E). From 
now on, we will use the notation 

A rh 2 = rh 2 — fh 2 c (33) 

to denote the deviation of the mass away from the thermodynamic critical point. 

In Fig. 7 we plot the full correlation functions £*-(z = 0) as functions of \i — j\ for various 
values of A?h 2 . In the insulating phase, the correlation functions decay exponentially as 
expected. The correlation lengtli tends to diverge as the critical point is approached, as 
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FIG. 7: The correlation functions t*- at the UV boundary (z = 0) as a function of' \i — j\ at 
different Aifi 2 for L = 13. (a) Logarithmic plot of at A m 2 > 0. The exponential decay becomes 
slower as the critical point is approached from the insulating phase. (b) t*j at A m 2 < 0 in the 
linear scale. In the superfluid phase, the correlation functions approach non-zero values in the large 
\i — j | limit, exhibiting the off-diagonal long-range order. 


is shown in Fig. 7 (a). In the superfluid phase, the correlation functions approach non- 
zero values in the large distance limit, exhibiting a long range order. As the critical point 
is approached from the superfluid side, the long range order disappears continuously as is 
shown in Fig. 7 (b). 

At the critical point, we expect that the correlation function decays in a power-law, 

■ 0(0)) ~ A. (34) 

with A = 1/2 for D — 3. In Fig. 8 (a), the correlation function is shown at A m 2 = 0. 
The curvature in the log-log plot indicates that the correlation function actually decays 
slower than an algebraic decay. This is due to the finite size effect which overestimates the 
correlations due to the periodic boundary condition. The finite size effect effectively pushes 
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FIG. 8: The UV boundary (z = 0) correlation functions t* 3 at Ara 2 = 0 for L = 5,7,9,11,13. 
(a) Log-log plot of the correlation functions. (b) Log-log plot of the correlation functions after 
constant pieces G(L) are subtracted out. (c) The exponent A (L) defined in Eq. (35) as a function 
of L at Ara 2 = 0 and Ara 2 = —3 (the line with Ara 2 = —3 will be discussed at the end of Sec. 
IV B. 2). In the thermodynamic limit the exponent approaches A(oo) ~ 0.55. (d) The offsets 
G(L ) defined in Eq. (35) as a function of L” 1 . 


the system at A m 2 = 0 to the superfluid side at finite L. As will be shown in Sec. IV B. 
3, the correlation function decays algebraically at a positive value of A m 2 that goes to zero 
only in the infinite L limit. To take this into account, we include a constant piece and fit 
the correlation function with 

t*j = G( L ) + |^ _ j|2A(L)> ( 35 ) 

as is shown in Fig. 8 (b). In the thermodynamic limit, the exponent A (L) approaches 0.55 
and G(L ) tends to vanish as is shown in Figs. 8 (c) and (d). We expect that the 10% error in 
the critical exponent can be made smaller by increasing the number of basis points that are 
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included to parameterize the radial profile of Uj{z) and t*j(z) in the nnmerical calculation. 
For details, see Appendix C. 


B. Bulk 

Now we set out to understand the different phases from the behaviors of the bulk fields 
in the (D + l)-dimensional space. 

1. Insulating phase 



FIG. 9: (a) Logarithmic plot of tjj as a function of z in the insulating phase at A rh 2 = 6 for 

L = 13. The curves from the top to the bottom represent the hopping fields with increasing | i — j\ 
from the nearest to the farthest neighbor hoppings. (b) Logarithmic plot of Uj(z ) as a function of 
|* — j | at different values of z. The fits with straight lines represent exponential decay. 

In Fig. 9, the scale dependent hopping fields are plotted in a deep insulating state at 
Arh 2 = 6 for L = 13. Although all further neighbor hoppings are generated in the bulk, 
they remain exponentially small in | i — j\ at all z. Furthermore, all hopping fields decay 
exponentially in z in the IR region. Due to the exponential decay both in |i —j\ and z, tij(z) 
is signihcant only for small \i — j\ and As the connectivity of space becomes weaker with 
increasing z, the correlation functions t*j(z) becorne smaller as is shown in Fig. 10. As 0 
increases, the overall magnitude of t*-(z) decreases exponentially while the rate of decay in 
\i — j\ is unchanged. In the large z limit, the space complctely loses its connectivity, and gets 
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FIG. 10: Logarithmic plot of t*j as a function of |i — j\ in the insulating phase at Am 2 = 6 
for L = 13. The overall magnitude decreases with z. The fits with straight lines represent an 
exponential decay with the same exponent for ali z. 
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FIG. 11: Logarithmic plot of the susceptibility x as a function of z in the insulating phase with 
Ara 2 = 1.5(a) and 0.5(b) for L = 13. Fragmentation scale is the crossover scale beyond which x( z ) 
becomes independent of z and becomes a constant whose value is determined by the on-site action 
Sq without hopping. For concreteness, we dehne zf to be the scale at which the second derivative 
of the logarithm of x( z ) is maximum. 


fragmented into isolated islands. In the RG language, the IR fixed point is described by the 
decoupled sites. One can quantify the scale of fragmentation in terms of the susceptibility 

x(z)=Yl t U z )’ ( 36 ) 

j 

which is proportional to ^ 2 in the insulating phase, where ^ is the correlation length. As 
z increases, the susceptibility decreases because the hopping amplitudes become smaller in 
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FIG. 12: Holographic phase diagrams for L = 13. The color contours represent the value of the 


susceptibility defined in Eq. (36). The darker the color is, the smaller the susceptibility is. zp is 
the fragmentation scale, and zh represents the horizon (see text in Secs. IV B 1 and 2). 

the bulk. In the large z limit, x saturates to the value determined by the on-site mass term. 
This is shown in Fig. 11. Although fragmentation is not a sharp transition, we can choose 
a convenient criterion for dehniteness: we define the fragmentation scale Zp to be the scale 
at which d is maximal. It represents a scale around which a crossover occurs from a 
connected lattice to fragmented space. In Fig. 11, zp is indicated in the insulating phase, 
which shows that zp increases with decreasing m 2 , as expected. Fragmentation occurs at 
longer distance scales when the systern is closer to the critical point. The background color in 
Fig. 12 represents the value of the susceptibility in the bulk for different values of A m 2 . The 
curve denoted by zf represents the fragmentation scale. The curve zp and other features in 
the superfluid phase (A rh 2 < 0) of Fig. 12 will be discussed in the following section. 

Near the critical point, the fragmentation scale Zp can be used to extract the critical 
exponent u, which dictates how the correlation length depends on the distance to the critical 
point, £ ~ (Am 2 )~ u . The susceptibility is given by x ~ £ 2 at the UV boundary. Inside the 
bulk, it falis off exponentially in z as xi z ) ~ £ 2 e“ 2z because of Eq. (30) before saturating to 
a fixed value of Aj. Therefore, the fragmentation scale is given by zp ~ ln£ ~ — v ln(Am 2 ). 
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FIG. 13: The correlation length £ plotted as a function of Am 2 on a log-log scale. The data is 
obtained from solving the UV boundary theory with lattice size L = 151. In the region where 
1 -C £ -C L/2, the scaling is close to the one expected in the thermodynamic limit, £ ~ (Am 2 ) , 

which is shown as a straight line. 

In Fig. 13 we show £ as a function of (Am 2 ), which shows v 1, the expected value for the 
interacting modcl we are studying in the large N limit. Therefore the fragmentation scale 
should behave as Zf ~ — ln(Am 2 ) in the thermodynamic limit. 

2. Superfluid phase 



FIG. 14: Hopping fields t^j in the superfluid phase at Am 2 = —3 for L = 13 plotted in the linear 
scale (a) and in the logarithmic. scale (b). The curves from the top to the bottom represent the 
hopping fields with increasing \i — j\ from the nearest to the farthest neighbor hoppings. 

In Fig. 14, we display the hopping fields inside the superfluid phase at Am 2 = —3 for 
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FIG. 15: (a) The log-log plot of the hopping fields as functions of |i — j | in the superfluid 

phase at A rh 2 = —3 for L = 13. From the UV boundary to a critical scale zh ~ 2.8, tij decays 
exponentially. The solid curves represent fits of the form in Eq. (37) with ijj(z) > 0. (b) The rate 
of the exponential decay ij)(z) from (a) plotted as a function of z. ip(z) continuously vanishes at a 
critical scale zh ~ 2.8. 


L = 13. Near the UV boundary, tij(z) decays exponentially in | i — j\ as is the case in the 
insulating phase. This has to be true even in the superfluid phase because only nearest 
neighbor hoppings are turned on at the UV boundary, and further neighbor hoppings are 
generated gradually in the bulk. The exponential decay of the hopping fields is well captured 
by 

-ip(z)\i-j\ 

‘«( 2 ) ~ y—yw < 37 > 

with 'ijj(z) >0 near the UV boundary as is shown in Fig. 15(a). However, the rate of the 
exponential decay, i]){z) becomes smaller as 2 increases, and vanishes at a critical scale zh- 
This is shown in Fig. 15(b). 

For z > zh, the hopping fields no longer decay exponentially but decay algebraically with 
a constant off-set, 

t ij{z) = C{z) + |——p^y, (38) 

where the constant off-set turns on continuously as 

C{z) ~ \z-z H \ pH (39) 

with an exponent /3h ~ 1 for z > zh- After the off-set is subtracted out the hopping fields 
exhibit pure power-law decays as is shown in Fig. 16(a). The z dependence of the off-set 
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FIG. 16: (a) The hopping fields Uj at Arh 2 = —3 after a constant off-set is subtracted out in the 

bulk with z > zh = 2.8 for L = 13. The straight lines represent pure power-law decays. (b) The 
off-set, C(z ) in Eq. (38) plotted as a function oi z — zh for L = 13. 


near zh is displayed in Fig. 16(b). 
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FIG. 17: (a) The off-set C(z) at A fh 2 = —3 for L = 7,9,11,13. (b) The log-log plot of the 
maximum value of C(z) as a function of L -1 at Arh 2 = —0.1,—1.5,—3. The lines represent fits, 
Cmax ~ L r,lm'z w ilb r(—0.1) ~ 2.7, r(—1.5) ~ 2.88 and r(—3) ~ 2.93, which suggests that C(z) 
vanishes at all Atti 2 in the thermodynamic limit. 


The scale Zh represents a ‘horizon’ at and beyond which the systern loses locality as sites 
are connected with each other through non-local hoppings. As a resnlt, the ‘coordinate 
distance’ |i — j | is no longer a good measnre of the ‘physical distance’ that dictates the 
actual correlations. Because of the non-local hopping, the physical distance between two 
sites effectively shrinks to zero. If we were studying a systern with a finite temperature with 
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a compact Euclidean time direction in such a way that its lcngth /3 would remain finite as 
L —> oo, we would first observe this loss of locality only in the Euclidean time direction. In 
that case, Zh corresponds to the usual horizon at which only the thermal circle shrinks to 
zero. I 11 the present case, the L —)■ oo limit corresponds to the zero temperature limit, where 
all directions are equivalent. As a resuit, the size of space shrinks to zero in all directions at 
the horizon. It will be interesting to understand what this generalized horizon corresponds 
to in Minkowski space. 

If the off-set C was indeed nonzero in the thermodynamic limit, it would irnply that 
any two sites remain coupled witli a non-vanishing hopping no matter how far they are. 
However, Fig. 17(a) shows that C systematically decreases as the system size increases. In 
order to understand the behavior of C in the thermodynamic limit, we need to do a finite 
size scaling analysis. Fig. 17(b) shows that the maximum value of C at fixed A m 2 decays 
as C max ~ -h with r = 2.81 ± 0.12 for —3 < A m 2 < —0.1 as the system size L increases. 
This implies that C goes to zero, and the hopping fields decay in a purely algebraic manner 
for z > Zh in the thermodynamic limit. 



Am^(L) 



L~ l 


FIG. 18: (a) The inverse of the horizon scale zh 1 as a function of L^ 1 with A m 2 = 

—3,—1.5,—0.1,0.1. The nonzero intercept in the large L limit implies that the horizon arises 
at a finite scale in the thermodynamic limit. (b) The critical mass determined from the location of 
horizon in the bulk, Am^(L), plotted as a function of L -1 . The line is a fit to an algebraic function 
of the form f(L ) = AL~ D with A ss 181 and D « 2.76. 


Although the off-set C vanishes in the thermodynamic limit, the location of the horizon 
Zh remains finite, as is shown in Fig. 18(a). In the deep superfluid phase, zh is more or 


32 






less independent of A rn 2 , but zh sharply increases as the critica! point is approached from 
the superfluid side. This is displayed in Fig. 12. For L = 13, the horizon scale diverges at 
A rh 2 = 0.2, and the horizon no longer exists for A m 2 > 0.2. From the location of the horizon 
in the bnlk, we define A rh 2 (L) to be the critica! mass for finite lattices, e.g., Am^(13) = 0.2. 
Although the critical point is at A m 2 = 0 in the thermodynamic linxit, Arh 2 (L) is different 
from zero for finite L due to the finite size effect. A finite size scaling in Fig. 18(b) suggests 
that A m 2 (L) indeed goes to zero in the thermodynamic limit. 
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FIG. 19: (a) k(z) from Eq. (38) as a function of z in the superfluid phase at Afri 2 = —3. It 

continuously decreases with increasing z inside the horizon, (b) L-dependence of k(z) at and inside 
the horizon at Ara’ = —3. 


In the absence of the off-set in the thermodynamic limit, what distinguishes the horizon 
at z = Zh from the region inside the horizon with z > Zh is the exponent k with which 
the hopping fields decay. Fig. 19(a) shows that n inside the horizon is smaller than kh 
defined at the horizon. The exponent k tends to approaeh no n zero values for all z in the 
thermodynamic limit as is shown in Fig. 19(b). 

Remarkably, the exponent kh at the horizon is independent of Am 2 as is shown in Fig. 
20. This suggests that the horizon is characterized by a universa! exponent associated with 
the decay of the hopping fields. The location where the decay of the hopping fields exhibits 
the universal exponent precisely coincides with the location where C turns on and ij} vanishes 
in finite lattices. It is interesting to note that the hopping fields exhibit universal ‘critica! 
behavior’ that is akin to a continuous phase transition. This shoulcl not be confused with the 
usual phase transition where the correlation length diverges as A m 2 is tuned from the insu- 
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FIG. 20: The log-log plot of the dynamical hopping fields tij as functions of \i — j\ at the horizon, 
zh , at different values of A rh 2 in the superfluid phase. The slope of the decay is the same at ali 
values of A rh 2 , which indicates the universality of k along zh- 

lating phase to the superfluid phase. The ‘transition’ that happens at z H is associated witli 
the divergence of the length scale for the renormalized hopping fields within the superfluid 
phase. 

As the ‘short-distance’ modes are integrated out, the effective action is renormalized by 
further neighbor hoppings. In the insulating phase, the further neighbor hopping fields 
remain exponentially small, and the renormalized action remains local at all scales. I 11 the 
superfluid phase, the System can not keep the locality as further neighbor hoppings are 
proliferated beyond the critica! scale zh- The presence/absence of the horizon can be used 
as a ‘holographic order parameter’ that distinguishes the superfluid/insulating phases. 

The non-locality in the bulk could have been avoided if one had allowed for spontaneous 
symmetry breaking by turning on an inhnitesimally small symmetry breaking field before 
taking the thermodynamic limit. In that description, one starts with a new vacuum with 
broken symmetry, and fluctuations around the new vacuum are described by the Goldstone 
modes. Here the spontaneous symmetry breaking is not allowed because the large volume 
limit is taken in the absence of the symmetry breaking field. The emergence of the non-local 
geometry is a sign that the System ca n n ot, maintain both locality and the full symmetry in 
the superfluid phase. 

As z iucreases further, t i3 (z) eventually decays exponentially in z while remaining flat 
in |i — j |. This lcads to fragmentation in the deep IR region even in the superfluid phase. 
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FIG. 21: The z-dependence of x in the superfluid phase at A fh? 


— 0.5(a) and — 1.5(b) for L = 13. 
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FIG. 22: L-dependence of zf in the insulating phase and the superfluid phase. The L axis is 
shown on a logarithmic scale. The non-zero slopes imply that zf diverges logarithmically in L 
in the superfluid phase. On the contrary, zf remains finite in the thermodynamic limit in the 
insulating phase. 

Fig. 21 shows that the susceptibility x becomes negligible beyond a fragmentation scale zp- 
The fragmentation scale Zf for L = 13 is displayed along witli the horizon scale zh in Fig. 
12. However, the fragmentation in the superfluid phase is a finite size effect unlike in the 
insulating phase. This can be seen from Fig. 22, where Zf increases logarithmically in L in 
the superfluid phase. This is in contrast to the insulating phase where Zp is independent 
of L. In the thermodynamic limit, zp — oo in the superfluid phase wliilc zp stays finite in 
the insulating phase. In this sense, the superfluid phase is characterized by the algebraically 
non-local geometry in the IR limit. 

In Fig. 23, the correlation functions t*j(z) are plotted in the bulk in the superfluid 
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FIG. 23: Log-log plots of t*- as a function of |z — j| at different z below and above the horizon, 
zh, in the superfluid phase at A m 2 = —3. 

phase. t*j(z) at each z is fitted to Eq. (35). The constant piece describes the long-range 
order, and the power-law decay originates from the Goldstone mode in the superfluid phase. 
The exponent A(z) remains constant as a function of z. In the thermodynamic limit, the 
exponent approaches A (L = oo) ~ 0.55 (Fig. 8 (c)). The correlation function in the bulk 
is insensitive to the divergence in the length scale of the hopping field across the horizon, 
which is at zh = 2.8 at A m 2 = —3. 



3. Critical point 
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FIG. 24: (a) A log-log plot of the hopping fields Uj as a function of |i — j | at several values of z 

at the critical mass Arh 2 = 0.2 for L = 13. (b) The logarithmic plot of as a function of z at 
A m 2 = 0.2. 
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2.5 3.0 3.5 4.0 4.5 5.0 5.5 Z 

FIG. 25: The logarithmic plot of 0(z) as a function of z for a larger lattice with L = 21 at the 
finite size critical point, A m 2 = 0.05. In the IR region, there is a ciear exponential decay in the 
form of 0(0 ~ e _1 ' 9z , which is shown as a straight line. 

In Fig. 24(a), we show the hopping fields Uj(z) as a function of | i — j | for L = 13 at the 
finite size critical point, A m 2 = Am 2 (13) = 0.2. The hopping fields are wcll fit by the form 
in Eq. (37). As z increases, the rate of exponential decay 0(0 decreases as is shown in Fig. 
24(b). Because of the invariance under the scale transformation at the critical point, 0( z ) 
is expected go to zero exponentially in the large z lirnit (it is noted that z is the logarithmic 
length scale). However, 0 in the IR region deviates frorn the expected exponential behavior 
due to finite size effect for small L (Fig. 24 (b)). To reduce the finite size effect, we pusli 
our system size to L = 21 to read off 0(0 at the critical point. As is shown in Fig. 25, 0 
indeed decays exponentially in the IR region. 

At the critical point, the hopping fields retain locality with the characteristic length scale 
1/0. In other words, the sites are more or less globally connected at length scales smaller 
than 1/0, and the theory is local only at larger length scales. However, we don’t expect 
there to be a sense of flat geometry at the length scale of 1/0 because the present theory is 
weakly coupled in the large N lirnit. 

As A/7j 2 is tuned across the horizon at a fixed z deep inside the bulk, tij shows the same 
critical behaviour as it does for increasing z with fixed A m 2 in the superfluid phase, which 
is discussed in Sec. IV. B. 2. Tliis is shown in Figs. 26(a) and 26(b). On the insulating 
side, the hopping fields decays exponentially, where the rate of the exponential decay (0) 
continuously vanishes as the critical point is approached. On the superfluid side, the hopping 
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FIG. 26: (a) The hopping fields tij at z = 4 plotted on a log-log plot as a function of |i — j\ for 

L = 13. Several values of Am 2 > Am 2 (13) in the insulating phase are shown, moving closer to the 
critical point line (Arh 2 (13) = 0.2). The solid lines represent fits with the form in Eq. (37), and 
the exponential component -0 is displayed for each curve, (b) The same log-log plot as (a), but at 
A m 2 < Am 2 (13). Here the constant off-set in Eq. (38) has been subtracted out to show that the 
remainder is a pure algebraic decay. 
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FIG. 27: C(z = 4) and ip(z = 4) for L = 13 from Eqs. (38) and (37), respectively, plotted as 
functions of A m 2 across the critical point inside the bulk. 

fields decay algebraically with the off-set C from Eq. (38) that continuously turns on across 
the horizon for Enite size lattices. The dependence of and C on A m 2 across the horizon 
at a fixed z is displayed in Fig. 27. Once the constant off-set is subtracted, the hopping 
fields decay algebraically inside the horizon. This is displayed in Fig. 26(b). Whilc the 
hopping fields decay with the universal exponent at the horizon, the exponent inside the 
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FIG. 28: k(z = 4) for L = 13 from Eq. (38) plotted as a function of Am 2 < 0 inside the bulk. 
horizon changes continuously, as is shown in Fig. 28. 
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FIG. 29: A log-log plot of i*- as a function of |f —j| at different z with fixed A?fi 2 = 0.2 for L = 13. 
The fields display an algebraic decay with the same power at all 2 . 

In Fig. 29, we show the correlation fnnctions t*j(z) at the finite size critica! point Afh 2 = 
0.2. They show purely algebraic decays of the form in Eq. (35) with G — 0. As is the 
case in the superfluid phase, A is independent of z and approaches A(oo) ~ 0.55 in the 
thermodynamic limit. 

4■ Holographic phase diagram 

Now we combine all the information to construet a holographic phase diagram in the 
bulk. In Fig. 12, the scale Zf marks the crossover beyond which the space is fragmented. In 
the insulating phase, Zp is largely independent of L : the fragmentation of space is a generic 
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feature of the insulating phase. On the other hand, zf diverges in the thermodynamic limit 
for A m 2 < 0. This is confirmed from the finite size scaling in Fig. 22. Therefore, there is 
no fragmentation outside the insulating phase in the thermodynamic limit. The superfluid 
phase is distinguished from the insulating phase by the presence of the horizon in the bulk. 
The horizon is characterized by the power-law decay of the hopping held, Uj(zH ) ~ \i-)\ K H 
with the universal exponent that approaches kh ~ 1 in the thermodynamic limit. At the 
critica! point, the hopping fields decay as Ujiz) ~ where ip{z), k(z) > 0 and 

i)(z) —* 0, k(z) —* kh as z —* oo. Fig. 1 summarizes the holographic phase diagram in the 
thermodynamic limit. 

V. SUMMARY AND DISCUSSION 

In summary, we derive the holographic action for the U(N) vector model regularized on 
a lattice. The bulk equations of motion are solved numerically for finite lattices in three 
dimensions. From the finite size scaling, we hnd that the insulating phase, the superfluid 
phase and the critical point exhibit distinet geometric features in the bulk. The IR geometry 
of the insulating phase is characterized by ultra-locality with a decoupled lattice. The su¬ 
perfluid phase exhibits a horizon at and beyond which the geometry becomes non-local. The 
critical point shows a local geometry with a characteristic lengtli scale that asymptotically 
diverges in the IR limit. 

Although the U(N) vector model is exactly solvable in the large N limit, the present 
holographic description allows one to understand the concrete model from a holographic 
perspective. This is a first step toward applying QRG to more non-trivial models whose 
Solutions at large N are not known (e.g. matrix models). An advantage of using finite 
lattices is that the finite size effect can be used in solving the problem using a simplihcd 
IR boundary condition numerically. We hnish by listing some open problems for future 
investigation. 

Analytic solution : Although the numerical Solutions provide a great deal of information 
on the behaviour of the bulk, it is desirable to have an analytic solution, especially at the 
critical point. One approach is to reduce the infinite set of equations of motion for the bi-local 
fields to a finite set by projecting the Solutions to the ones constrained by the numerical 
solution. In general, it will be of interest to better understand theoretical structures of 
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bi-local (or multi-local) field theories. 

Finite temperature : The Euclidean theory in D dimensions can be viewed as the (D — 1)- 
dimensional quantum theory at zero temperature in the imaginary time formalism. One can 
turn on hnite temperatures by making the size of the thermal circle finite. In tliis case, the 
hopping helds in the temporal direction will be different from those in the spatial directions. 
At the horizon, the non-locality is expected to develop only in the temporal direction but 
not in the spatial directions. In other words, the size of the thermal circle will shrink to zero 
while the spatial area remains nonzero at the horizon. 

Application to Ferini surfaces : Having understood the insulator to superfluid phase 
transition in the bosonic model, one can try to study the fermionic counterpart. In the 
fermionic System, the superfluid phase is replaced by an itinerant state with a Fermi surface. 
It will be of great interest to understand how the Fermi surface manifests itself in the bulk 
geometry [34-38] . 

Full (D + l)-dimensional diffeomorphism invariance in the continuum limit : As is dis- 
cussed in Sec. II, the holographic action for the lattice model possesses a subset of the 
full diffeomorphism invariance of the continuum space. At each O-dimensional slice with a 
constant z, the local permutation symmetry is present. Given that the lattice provides a UV 
complete theory which flows to the conformal field theory in the continuum (long-distance) 
limit, it is expected that the holographic theory for the lattice model recovers the continuum 
theory with the full diffeomorphism invariance in the IR region. However, this needs to be 
understood more explicitly. 

Connection to the liigher spin theory : In the continuum limit, it may be possiblc to relate 
the equations of motion in the bulk with those of the liigher spin gauge theories proposed as 
the holographic dual for the vector models[33, 39-44]. However, the connection is not ciear a 
priori because the form of the bulk theory is sensitive to the specikc regularization scheme. In 
particular, we believe that the 1 /N corrections cannot be included in the present formalism 
of the liigher spin theory where the quartic interaction is implemented only through the UV 
boundary condition[33]. This is because the quadratic theory with fluctuating sources is 
ill-dehned unless one keeps the double-trace operator in the bulk, as is discussed in Sec. II. 
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Appendix A: Derivation of Eq. (14) 


Here we show the intermediate steps in going from Eq. (13) to Eq. (14). We can rearrange 
§( 0 )" 

in the following way: 


S ( °)" = NS UV d 0) l + S L [0, 0*] + S H 


(j), 0* 


+ §i 


0 , 0 *, 0 , 0 * 


A 




where 

§ L = m 2 (<£* ' 4>i) + 

i i ij 

§ h = Y1 Vi (<t>* ■ 0 *) , 
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S/--E 0 >e “ [ 0 ? ■ <0 + (0* • 0,) + (0* ■ 0, 


(0* ■ 0,) 2 - v ( “ s ‘ >+ Y)‘ k (0* ■ 0,). 
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N ^ 


dz 


0?-0i + #■& + #> 


+2 (0* • 00 | (0? ■ 0;J + [ 4 >* ■ 0ij + (^0* • 0j 

+2 ^0* ■ 0;) | ^0* ■ 0i) + ^0* ' 0i) } + 2 (0* ■ 0;) ^0* • 0i) 

Sl [0, 0*] is the bare action for the low energy fields 0 and 0*, S// 0, 0* is the bare action 
for the high energy fields 0 and 0*, and S/ gives the mixing between the low and the liigh 
energy fields. Since a■ 1 ^ and m are both 0(1) in dz, if we rnake dz infinitesimal we have 

1 e 2a ? >dz - 1 2 a^dz 


P-i 


O(dz). 


rn m- 

Therefore, to derive the bulk action that is continnous in the holographic direction we need 
to take into acconnt only ternis that are linear in 1 /jl 2 . Keeping this in rnind we integrate 
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out the high energy fields <fi, (f>* using only their bare action. For this we write Z as 

Zoc f D(/)D0*D^ (o) D^* (o) e■ s l O),, e~ As l 0,,, , 


where 


sf y ' = m uv [(S">,d 0) l +s t [<A,0 


and 


where 


e-^i 
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= / D(t)V(j)* e~ $H e~ $I = yK(e- Sl )s H , 


-s i 


= / D(pT)(p*e 


* S f 


coV' 

is the normalization factor. We need to calculate the correction to the action AS) ; to first 
order in dz. To make the calculation more tractable we rewrite 

=^(l + f; hf((S,)')i, 

v U p - 


and renarne the ternis from the sum 


f-lfp 

^ = L -r((S/) p )* 

p\ 


The goal is to calculate every to order dz , then re-exponentiate to get AS) J to order 


dz. is given by 


2?x = 


2 dz 


m z 




(Al) 


The first term of Eq. (Al) comes from tracing out the high energy modes from the hopping 
term, while the second term comes from the renormalization of the low energy mass due to 
integration of the high-low energy 4-boson interaction. The two contributions are illustrated 
in Figure 30. ^ is given by 




1 2 dz \ > ( 
2m^ ^ ai 


(i) 


2 + 


8A 2 


jk 

4A 


N 2 


77 (*y 0) (tf ' • &)) (0* • 0*) 


(A2) 


45 












/(0) A 

L ii N 


(a) 


(b) 


FIG. 30: The diagrams contributing to at order dz. (a) and (b) represent the first and second 
terms of Eq. (AI), respectively. The dashed propagator indicates the high energy modes while the 
solid propagator indicates the low energy modes. 


The three contributions are illustrated in Fig 31. The first one represents fusion of hopping 
links to generate further hopping links (a). The second represents the fusion of two 4-boson 
vertices to generate a 6-boson vertex (b), while the final two terms refer to the fusion of 
the 4-vertex witli the hopping term (c). There are no more terms in any that are linear 
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FIG. 31: The diagrams contributing to ^2 at order dz. (a), (b), (c) represent the first, second and 
third terms of Eq. (A2), respectively. Dashed and solid propagators have the sarne meaning as in 
Fig. 30. 
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in dz, which implies that AS^ 0) = — (3?\ + 3f 2 ) + 0(dz 2 ). This finishes our calculation of 
§i°^ + AS^ , which in Eq. (14) we call Ih 0 )". 


Appendix B: Large N solution for U(N ) vector model on a lattice 


Here we review the solution of the U(N ) vector lattice model at N —y oo [45]. Our initial 
partition function at any given scale is 



We introduce a Hubbard-Stratanovich field a in the following way 



Redehning cr* = —fer* gives 



We want to replace <7* with its saddle point value at N —> oo, which is assumed to be a 
site-independent constant <7. To find cr we integrate out the and <fi* fields after going to 
momentum space 


J d^dd) e~y^- tk+m2+2<TMk \ 2 } = fjj 


TtV \ Ntrlogf nV -'i 

__ — e y V m 


m 2 — ffc + 2a 


The effective action for a is now (ignoring constant ternis) 

Seff = — V—a 2 + N tr log (in 2 — 4 + 2cr) , 
A 


and the saddle point equation becomes 

dS eff 


da 


-2 V^a + Nj2 2 


m 2 ~ tk + 2cr 


= 0 


_ _ 1 ^ A 

T/ rr > 2 — 


V ^ m 2 — tk + 2cr ’ 


which is the self-consistency condition for the rnean field value <7, Eq. (26). 
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Appendix C: Numerical Method 


Here we discuss in detail the method we use to solve the equations of motion (EOM) 
of Eq. (21). The equations themselves are first order and seemingly simple, however each 
equation for tij and £*■ involves a sum over ali other fields t qp and/or t* p . The number of 
EOM grows as the volume V, and the number of terms in the sum also grows as V, leading 
to a very large number of terms. The equations are also non-linear. We use a tau spectra! 
collocation method [46-48] combined with a Newton-Raphson algorithm [49] to solve this 
System. 

Each bilocal field is decomposed into a truncated polynomial series 

n n 

tii(z) = Y J <T a (x(z)), f ij (z) = J2^T a (x(z)), (Cl) 

a =0 Q !—0 

where T a (x(z)) is the a-th Chebyshev polynomial of the first kind. Here x(z) rnust live on 
[—1,1], which is the domain on which the Chebyshev polynomials are orthonormal. The 
domain of the RG scale z is [0, z*], which implies that x(z) = — 1. The number n is called 

the order of the tau spectral approximation, and n+1 is also the number of collocation points, 
i.e. the number of points on the x/z- axis at which the spectral Chebyshev representation 
agrees exactly with the original function. n is also a measure of the accuracy of the solution; 
for smooth functions, one expects the error in the solution to decrease exponentially with 
n. In our case we take the collocation points to be the Chebyshev Gauss-Lobatto (CGL) 
points x a = cos(^). For each lattice size we take n — 20. 

We solve our equations for the coefficients a,b instead of the full functions t,t*. This 
allows us to solve one large (non-differential) matrix equation rather than rnany differential 
equations. For this we rewrite the EOM and boundary conditions in terms of the new “vec- 
tors” a t j and bij, where the {ij} labcl the lattice link and the vector labcl is the Chebyshev 
label (labeled as “a” in Eq.(Cl)). 

The rewritten EOM have 4 types of terms. The first is a constant, which is represented 
as a vector with that constant as the first component (a = 0) and all other components 
zero. The second type is a lincar term, which is just CijUj (or Ci 3 t*-) and it becomes a\j 
(or b^). The tliird type is the differential operator. We denote it by the matrix L and 
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for the CGL points it is given by 


if /3 > a and /3 + a = odd 
otherwise 


-^a/3 


{ 2 2/3 

z* c(a+l) 

0 


where c(y) = 1 + <5 7i0 . The last type of term is the quadratic term, e.g. t pq . We choose to 
express it in spectra! form via a procednre described in Sec. 3.4.4 of [47]. The first step of 
the procedure is to expand the sum to order m = 2 n, leaving ali the coefficients zero that 
are higher order than n, i.e. a ^ = a\ = 0 for n < a < 2 n. Then we perform an inverse 
discrete Chebyshev transform to obtain the real space coefficients for every field, which are 
given by 


~Z7 

a p = 


E 

CK—0 


a l Jlcos 


7rf3a 


m 


We multiply the real space coefficients to obtain the coefficients of the product in 
s‘p Pq = a l p ■ a P p, for /3 G {0,m}. We now perform an inverse discrete Chebyshev 
back to spectra! space, which is given by 


real space 
transform 


where 


Finally, 


Jj,pq 

d a 


m 


E 


s /3 


2 

m c a C /3 


cos 



> 


j 2, 7 = 0, m 

I 1 7 G { 1 , m — 1 } 


tij(z)tpq(z) = S a ,P9 T a(x(z)) 

a=0 


is the spectra! representation of the product of tij(z) t pq (z). 

Now it is straightforward to express the EOM and boundary conditions in spectral form 
as a set of algebraic equations. The number of equations is equa! to the total number of 
ficlds tij and £*■ times n+ 1 (the number of collocation points), which we enumerate and put 
into vector form R(aij, b iq ) = 0. Finally, we combine ali of the {<%, bij} into one big vector c 



(C2) 
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As R is quadratic in c, we cannot solve it directly via a matrix inversion. We employ an 
approximation method, namely a Newton-Raphson method, to solve Eq. (C2). For this we 
linearize the equations 

= 0 , 

CO 

and solve for c). This process is iterated w times to arrive at a solution c w , which satisfies 

L 2 = — cVi)a < e - We choose e = 10 —11 , which is small enough to indicate that 

V OL 

the Newton-Raphson method has converged up to floating point round-off. 

An important question is that of convergence: when does the Newton-Raphson method 
converge, depending on the initial starting vector cq, and how many iterations does it take? 
The convergence domain that we hnd is summarized below (the reader should keep in rnind 
that in tuning rh 2 = m 2 — tu we change t lt and keep m 2 constant (in this case m 2 = 25)). 

The solution in the deep insulating phase is available from the perturbative analytical 
solution of Eq. (29). From this starting point one can obtain the solution at a certain 
distance away from the critical point on the insulating side, say A m 2 = 10. From this 
solution the Newton-Raphson method converges quickly (within 10 iterations) for any point 
A/ 7 i 2 > 0. Also, it converges at A m 2 < 0, however more slowly (« within 20 iterations). 
This gives us the Solutions at and just beyond the critical point on the superfluid side, which 
are of the most interest physically. Starting from a solution at 10 3> A m 2 > 0, we see a 
faster convergence to these desired points. Accessing points in the deep superfluid state 
however is more difficult. To dehnitively converge to a point A m\ <C 0 one has to start witli 
A/ 7 /q = Am 2 + Stu, where 8t tl depends on the lattice size L. For L = 13 we find that St tl m 1 
is a good choice. The slower convergence in the superfluid phase is due to the sensitivity of 
R(c) in Eq. (C2) on c. 

In summary, the numerical code can be parallclized when hnding Solutions on the in¬ 
sulating side, as well as at and near the critical point on the superfuild side, but hnding 
Solutions in the deeper superfluid phase needs to be done in serial. 

Appendix D: Comparison between the general and on-site boundary conditions 

In this appendix we check that the on-site boundary condition (BC) gives the same 
Solutions to the EOM as the general BC. The Solutions for the L = 3 lattice with the on-site 


R(co) + (Ci — Cq) 


dR 

dc 
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BC are plotted in Figs. 32 and 33. We can see that the fields look identical to those of Figs. 
4 and 5, and after comparison we find that they indeed lic on top of each other. 


*(0,0,0) ^(0,0,1) 



f (0,l,l) *(U,1) 



(c) (d) 

FIG. 32: The same plots as in Fig. 4 with the difference being that the IR BC used is the on-site 
one of Eq. (31), instead of the general BC of Eq. (27). 


However, besides the bulk shape it is also important to compare the exact IR values to 
those obtained using the general BC, especially when the theory starts out in the superfluid 
phase at the UV. Tables I and II summarize this comparison for two values of frt 2 = 25 
and —15, corresponding to the deep insulating and deep superfluid phases, respectively. In 
both phases the t*^j(z*) values are the only ones that differ at all between the two BCs. For 
the on-site BC they are zero always, which is what we set them to (10 -15 is the machine 
precision), and for the general BC this is not the case. However, for the general BC 

are stili extremely small when compared to t^z*) in both phases, which is why it works to 
approximate them by zero. This explains the validity of this approximation. 
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(c) (d) 

FIG. 33: The same plots as in Fig. 5 with the difference being that the IR BC used is the on-site 
one of Eq. (31), instead of the general BC of Eq. (27). 


m 2 = 25 

tu{z*) 

)| 

Uz*) 


general BC 

0.08 

< 1.3 x IO" 7 

0.04 

< 2 x 10" 10 

on-site BC 

0.08 

< 1.3 x IO" 7 

0.04 

< IO" 19 


TABLE I: Comparison at z = z* of general and on-site BCs for fh? = 25 (deep insulating phase). 

In Fig. 34, we show the hopping fields t l3 {z) for larger lattices L = 5, 7, 9 obtained using 
the on-site BC. They look qualitatively similar to those in Fig. 32. 
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m 2 = —15 

tu(z*) 

\tijLj(z )| 



general BC 

0.079812 

w 0.00071 

0.0399997 

« 1.14 x 10“ 6 

on-site BC 

0.079812 

w 0.00071 

0.0399997 

< io - 15 


TABLE II: Comparison at z = z* of general and on-site BCs for rh 2 = —15 (deep superfluid phase). 

Appendix E: Value of m 2 in the L —> oo limit. 

The discrete version of our action in Eq. (2) has a critica! point in the thermodynamic 
limit, m 2 , which can be computed exactly in the large N limit. Using Eq. (26) a,t z = 0 
(in the thermodynamic limit the sum becomes and integral), we compute m 2 for the values 
of tij and Jijp q which we set at the UV boundary: Uj is zero beyond nearest neighbors and 
Jijpq = 0. At z = 0 the hopping fields in momentum space are 

= 'Yh = *(o,o,o) + 2 tnn ( cos(k x ) + cos(k y ) + cos(k z )), 

( l,m,n ) 

where t nn = 1 are the nearest neighbor hoppings. At the critical point the denominator of 
the integrand in Eq. (26) is zero exactly for k — 0, since t% = Q > i We can plug in the 
value of cr c for which this is true on both sides of Eq. (26), which gives 

1 (f 2 n _ f 2n f 2n f 2n _A_ 

2 ( ^(o,o,o) m c + ) y o J o J o ( 27 t) 3 6 — 2 (cos(k x ) + cos(k y ) + cos(k z ))' 

Reverting back to the notation of t(o,o,o) = tu, we obtain 

fh 2 c = (m 2 — tu) c = 5.49454. 
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*( 2 , 2 , 2 ) 


khiD 



(a) 


(b) 


^(2,2,2) ^(3,3,3) 



(c) (d) 

f (3,3,3) ^(4,4,4) 



(e) (f) 

FIG. 34: Same plots as Fig. 32 for larger lattices L = 5 (a,b), L = 7 (c,d), L = 9 (e,f). In all 
plots, rh 2 are chosen to be 15, 10, 7.5, 6.5, 6, 5.49454, 5, 4.5, 4, 3.5, 2.5, 0, —5, and the arrow frorn 
Fig. 32 applies here as well. The critical point (rh 2 = 5.49454) is again denoted by think curves. 
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